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Measurement of neutral pions and direct photons are closely connected experimentally, on the 
other hand they probe quite different aspects of relativistic heavy ion collisions. In this short 
review of the n'^ results from the PHENIX experiment at RHIC our focus is on the -integrated 
nuclear modification factor, its energy and system size dependence, and the impact of these re- 
sults on parton energy loss models. We also discuss the current status of high pr and thermal 
direct photon measurements both in p+p and AuH-Au collisions. Recognizing the advantages of 
measuring not only the "signal", but also all the "references" needed for proper interpretation in 
the same experiments (with same or similar systematics) we argue that RHIC should regularly 
include d+A and even d+d collisions into its system size and energy scan. 
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Figure 1: "Historic" nuclear modification factors Raa vs pr in 200GeV/A Au+Au collisions. Squares: 
direct photons (no apparent suppression Q, published 2005). Triangles and circles: preUminary results on 
n° and rj, also from 2005. 

1. Introduction 



One of the predicted signals of the formation of a quark-gluon plasma (QGP) in relativis- 
tic heavy ion collisions was the in-medium energy loss of hard scattered partons, leading to "jet 
quenching" and manifesting itself in a drastic reduction of high pT hadron production in relativis- 
tic heavy ion collisions, as compared to the production rate in p+p scaled by the nuclear thickness 
function Taa to account for the increased probability of hard scattering. The effect has been char- 
acterized by the nuclear modification factor Raa defined as 
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for an arbitrary particle type a. In the absence of nuclear effects Raa = 1, otherwise it is 
enhanced {Raa > 1> hke in the "Cronin-effect") or suppressed (Raa < 1 like in jet quenching.) 
The first major discovery at RHIC was that in central Au-i-Au collisions hadrons were strongly 
suppressed [jl], ^ (see Fig. [T]), in accordance with the predicted energy loss of partons in the hot, 
dense medium. Lack of suppression in d-i-Au collisions [^] (where we assume no such medium is 
formed) supported the conclusion that the jet suppression was indeed a final state effect, as opposed 
to some drastic changes in the initial state (parton distribution functions). The first observation 
that direct photons were not suppressed, not even in the most central Au-i-Au collisions |Q] made 
this first, qualitative picture coherent, by showing that the crucial step in forming Raa - scaling 
p+p cross sections by (7^) to obtain the expected Au-i-Au yields - makes sense: photons will not 
interact with the medium (a^m << GCs), so if the {Taa) scaling is the right concept, the direct photon 
Raa should be unity - and indeed, it is (see Fig |I|). 

While very useful and providing many early insights, causing a flurry of theoretical activities 
rarely seen before, the above picture turned out to be only rudimentary. As RHIC experiments 
moved "from the discovery phase to the exploration phase" more precise data became available. 
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with higher statistics, extended pr ranges and more complex physics quantities (like two- and 
three-particle correlations, azimuthal anisotropics, etc. Equally important, RHIC launched a sys- 
tematic energy and colliding species scan. This made measuring excitation functions in the same 
collider and within the same experiments possible. The collider environment ensures that multi- 
plicities at midrapidity rise only slowly with ^/s^n, and measuring in the same experiment means 
that most systematic errors will be similar - even if individual measurements are not very precise, 
the evolution with energy and/or system size can be established quite accurately. Increased accu- 
racy of the data made comparisons with theories more meaningful, put some constraints on free 
parameters [Q] and gave the first experimental indications at what region of the phase diagram the 
transition to the strongly coupled quark gluon plasma might occur. Almost needless to say that in 
the process our initial, relatively simple picture became richer (much more complicated). 

2. Pion suppression systematics 

So far four aspects of the ^-integrated Raa systematics have been studied at RHIC. First, the 
magnitude and shape of Raa{pt) at the highest available transverse momenta - namely, whether it 
remains constant, consistent with a constant fractional energy loss picture, or rises with increasing 
Pt, as predicted by most parton energy loss models These studies gave the first quantitative 
constraints on free parameters of models, like initial gluon density or the transport coefficient of 
the PQM ^ model = /i^/A) = 13.2+3 2GeV^/fm, where jj. is the average transverse momentum 
transfer to the medium per mean free path A. PQM - which includes radiative loss, a static medium, 
no initial state multiple scattering and unmodified PDFs - predicts a slow rise of Raa with pj, as 
do other models with very different assumptions In light of this it is interesting, that the data 
are best described by a constant fit which in turn for the power law spectra of high pj pions would 
imply a constant fractional energy loss (constant Stoss = ^Pt / Pt fractional shift of the spectrum). 
These statements are based on data taken in RHIC Run-4 (2004), which still have large errors on 
the high px points, results from the much higher statistics 2007 data should be available soon. 

Second, the evolution of Raa with collision centrality has been studied. One very interesting 
and non-trivial observation was that the shape of the high px Tl''' spectra is virtually unchanged 
from p+p to the most central Au-i-Au collisions. PHENIX found that the exponent of the power 
law function fitting the high pj part of the spectra (<=<: p'^) varies from n = —8.22 it 0.09 in p+p 
to n = —8.00 ±0.012 in the most central (0-5%) Au-i-Au collisions. This combined with the fact 
that Raa in the most central collisions is essentially flat (constant) means that Raa is constant at 
all centralities to the best of our current knowledge which in turn means that Raa integrated above 
a certain pr is a useful quantity to characterize the onset and evolution of jet quenching. This is 
shown on Fig. Fitting the points with a function motivated by the constant fractional energy 
loss {i.e the energy loss increases with pt) Raa = (1 — ^0^partT~'^ the exponent a is found to be 
a = 0.58 ±0.07 for pT > 5GeV/c and a = 0.56 ±0.10 for pr > lOGeV/c which is in reasonable 
(albeit not perfect) agreement with the a ^2/3 predictions of the PQM [^] and GLV ^ models. 
Note that in the "constant fractional energy loss" picture Raa depends on the slope (exponent n) of 
the power-law pt spectrum which in turn depends on the per-nucleon collision energy ^/sNN. 



Therefore, the third issue is what can be learned from the energy dependence of Raa, and the 
first answer is given on Fig. |3| which shows the nuclear modification factors measured at RHIC in 
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Figure 2: Integrated Raa for n'^ as a function of collision centrality expressed in terms of Npan- The last two 
points correspond to overlapping centrality bins 0-10% and 0-5%. See text for an explanation of the fits. 



the most central Cu-i-Cu collisions at ^snn =22.4,62.4 and 200GeV. Instead of suppression the low- 
est energy data show a substantial enhancement of Raa in hne with the well-known "Cronin-effect" 
(and within errors this enhancement is independent of centrality [^) whereas at y/s^N =62.4GeV 
the pions are already significantly suppressed. Of course it does not follow that the onset of jet 
quenching - presumably due to QGP formation - happens necessarily between y^S^v =22.4 and 
62.4GeV; all we can say is that the quenching mechanism overtakes the Cronin-effect in this re- 
gion. Being a consequence of the initial smearing of partons the Cronin-effect is present already 
in cold nuclei and could be unfolded from the yS^-dependence of Raa if P+A {d+A) data were 
available from RHIC in the entire available collision energy range. In the author's opinion taking 
d+K data (along with A-i-A and p+p for reference) should be part of the planned "energy scan" 
at RHIC. - Since pions are strongly suppressed in Cu-i-Cu already at ^/s^n =62.4GeV, it is not 
surprising that they are also suppressed in Au-i-Au, in fact, almost as strongly as at 130 or 200GeV. 

Finally, the fourth issue is how does Raa depend on the (size of the) colliding systems at any 
given energy. A comparison of ^/sMN =200GeV central Cu-i-Cu colhsions to mid-central Au-i-Au 
(such that the number of participating nucleons Npan is the same in both cases albeit the geometry 
is not) shows similar values of Raa- It would be interesting to know which is the lightest, symmetric 
heavy ion system where jet quenching at ^/sj^ =200GeV would disappear. 

So far we discussed only the ^-integrated Raa- For any given dataset neglecting azimuthal 
dependence results in higher precision of the measurement itself (better statistics and no system- 
atic errors related to the determination of the reaction plane) but it also hides any information 
encoded in azimuthal anisotropics of Raa due to the overlap geometry of the colliding nuclei. Such 
anisotropics clearly exist: after all, some of the most spectacular observations at RHIC are the large 
"elliptic flow" of hadrons and the quark number scaling thereof. Indeed, when studying /?^(A<I>) 
where AO is the azimuthal angle with respect to the reaction plane, determined event-by-event, we 
observe a rich structure [ IC]: jet quenching varies substantially with the average pathlength of the 
parton in the medium. Back-to-back dihadron (or more ambitiously, dijet) correlations, and par- 
ticularly the "golden channel", photon-hadron (photon-jet) correlations have also been advocated - 
and now measured - as tools of "jet tomography". It has been argued | [TT| ] that single hadron spec- 
tra at high pj are dominated by surface emission, therefore, Raa is barely sensitive to the medium 



4 



Pion and photon production in heavy ion collisions 



Gabor David 



< 
< 



- Vitev,\js^ = 2: . 

- no energy loss 



■ Ws^ = 22.4 GeV 



CC 




0.5 



1.5 



2 





Cu+Cu, 0-10% most central 



■ Vitev, 22.4 GeV, 130 < dN^/dy < 185 

□ Vitev, 62.4 GeV, 175 < dN^/dy < 255 

□ Vitev, 200 GeV, 255 < dN^/dy < 370 



■ o • 



ilf 







p^^(GeV/c) 



Figure 3: Nuclear modification factors for the most central Cu+Cu collisions [^] for ^snn =22.4,62.4 and 
200GeV compared to theoretical calculations [|l2[ . 



deep inside the collision volume. On the other hand back-to-back correlations force one of the 
partons to go through the medium (with the exception of "tangential" emissions) thus the dihadron 
suppression Iaa is more sensitive to the input parameters of theories like initial energy or gluon 
density. This is certainly true as long as the uncertainties of both theory and experimental data 
are negligible: a multidifferential quantity always provides more constraints. However, theories 
usually have uncertainties (even if often unstated, regrettably), and experimental data always have 
them and it is hard to imagine a situation in which using a given dataset a multidifferential quantity 
could be measured more precisely than a single differential. Just the opposite. The quantitative 
question then is: does the increased sensitivity on the theory side compensate for the inevitable 
loss of precision in the data? The answer varies from case to case. 

3. Direct photons 

Direct photons - defined as those not originating from final state hadron decays - offer unique 
possibilities to study heavy ion collisions. First, there are physics mechanisms that produce di- 
rect photons at each stage of the collision, thus they have the potential to provide information on 
the entire evolution, "history" of the collision. Equally important, once produced, they leave the 
interaction region mostly unaltered (a^,,, << a,), so the information they carry is unbiased. Unfor- 
tunately, the very same properties make them challenging to measure (photons from hadron decays 
provide a large background, particularly at low pj) and their "message" is hard to decypher be- 
cause it is difficult (if not impossible) to deconvolute the contributions of individual mechanisms 
to the observed direct photon spectrum. 

On Fig. [T] we already have seen the first measurement of direct photon Raa in central Au-i-Au 
collisions at ^snn =200GeV, showing that in the region where photons from primordial hard scat- 
tering were expected to dominate {pj >5-6GeV/c) Raa is indeed consistent with unity. This was 
consistent with the picture that both the PDFs and the production cross section are unchanged, the 
only difference between p+p and Au-i-Au as far as hard scattering is concerned is the increased 
probability of an otherwise rare process; this increased probability is properly described by a 
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Figure 4: Direct photon Raa in central Au+Au collisions. The denominator is a fit to the PHENIX pre- 
liminary Run-5 p+p data. Theoretical curves are LO calculations from [[l5[|. The solid line shows the pure 
isospin effect, the dash-dotted line is the combined effect of isospin and (anti)shadowing (EKS), finally the 
band at the botoom combines isospin and antishadowing with photon quenching. 



straightforward geometric overlap integral (the nuclear thickness function Taa), and the photons 
produced escape the collision region unaltered even if a dense, strongly interacting medium is 
formed which quenches parton jets. Although we have higher quality data now (in an extended 
range) and our understanding of the processes generating photons is much more nuanced, we 
should emphasize that at a very basic level, in first approximation, qualitatively the above picture 
is still correct. 

What are the currently known/assumed sources of direct photons in p+p and Au-i-Au colli- 
sions? At leading order two processes generate photons, the quark-quark annihilation (qq —>■ gj), 
which is suppressed at high pr due to the lack of valence antiquarks, and the dominant quark-gluon 
Compton-scattering {qg qy). Note that in the latter the photon balances the momentum of the 
original quark (modulo the initial kj smearing of the partons) and it comes out freely from the 
medium, whereas the quark will lose a significant fraction of its energy. Therefore, in back-to-back 
photon-jet correlations the photon "calibrates" the energy scale of the jet, which is important in 
"jet tomography" as well as in studying possible modifications of the fragmentation function in the 
presence of a medium. The probe itself is very clean and well understood - unfortunately the rates 
are quite low, due to the ttg^ coupling. 

At higher order fragmentation (or Bremsstrahlung) photons contribute both in p+p and heavy 
ion collisions. In principle Compton and fragmentation photons are distinguishable (at least in 
p+p) because fragmentation photons should have enhanced, jet-like activity in a small cone around 
them (tracks or energy deposit that exceeds the level justified by the underlying event), whereas LO 
Compton photons should have none in their neighborhood: they are "isolated". There are predic- 
tions that at high pj fragmentation photons will be a larger fraction of all direct photons in A-i-A 
collisions than in p+p. This is another reason why is it so important to understand p+p processes as 
well as possible: they will serve as reference when we try to find new phenomena (new sources of 
direct photons) in A-i-A. First direct photon cross-sections in =200GeV p+p collisions at RHIC 
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Figure 5: The fraction of the direct photon component as a function of pr in (a) p+p and (b) Au+Au (min. 
bias). The error bars and the error band represent the statistical and systematic uncertainties, respectively. 
The curves are from a NLO pQCD calculation. 



were published in [|Tj] followed by [14] with extended pr range, small systematic errors and de- 
tailed comparison to NLO pQCD calculations. Above pr > 5GeV/c where the systematic errors 
of theory and experiment are comparable, the data are well described by the theory, including the 
fraction of isolated photons 

Is p+p the right reference at all? For pions it certainly is due to the isospin symmetry of 
the protons and neutrons of the nucleus. The same is not true for direct photons: due to the 
coupling the photon cross section is proportional to the sum of the squared quark charges Le^ 
which is different for p and n. This "isospin effect" causes a trivial shift of the direct photon 



Raa [115[], depending on the nucleus and to some small extent on centrality; in the absence of any 
other medium or nuclear effect Raa should go about 10-15% below unity in Au-i-Au collisions (see 
the solid curve on Fig. ^ which shows the latest, still preliminary PHENIX results for photon Raa 
in the most central Au-i-Au collisions). Note that this trivial effect may either mask or enhance other 
mechanisms affecting Raa- The clean way to generate direct photon Raa would be to compare to 
a proper mixture of photon cross sections in p+p, p+n and n+n collisions, which could be readily 
available at RHIC from d+d collisions, since these three types of collisions could easily be tagged. 

A recent overview of the possible (additional) photon production mechanisms in heavy ion col- 
lisions is given in [|l6|]. At medium pj (corresponding to the 0.1 < x < 0.2 region) anti-shadowing 
should enhance the photon yield. The large parton energy loss (deduced from tt" suppression) 
means that fragmentation photons could be suppressed in the entire known pj range, unless the 
fragmentation functions themselves change in the presence of a medium (which is hard, but not 
impossible to measure in Au-i-Au). On the other hand the jet-photon conversion mechanism may 
increase high pj photon production. In this process a hard scattered quark interacts with a gluon 
or antiquark of the medium and the outcoming photon carries essentially the full momentum of 
the original parton (hence the name). Disentangling the individual contributions of these processes 
will not be easy, although we should point out that some mechanisms produce isolated, others give 
non-isolated photons, and their azimuthal asymmetries are also different (for instance jet-photon 
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conversion should produce isolated photons with a negative "elliptic flow" V2 < which is quite 
unique since all other particles tested so far at RHIC exhibited V2 > 0). 

Low pt real photons (thermal region) are very hard to measure due to the overwhelming 
hadron decay background, however, "quasi-real" photons are accessible as low mass e^e^ pairs. 
Using this technique PHENIX measured the fraction of direct photons in the total inclusive photon 
spectrum [|l7|] (see Fig. ||) and found that while in p+p the results are consistent with NLO pQCD 
calculations (there are no unexpected sources), in Au+Au there is a statistically significant 10% 
excess at low pj. 

In summary, after the early, exciting discoveries RHIC experiments entered the era of explo- 
ration. Data became more precise enabling detailed comparisons to theories and setting meaningful 
constraints on their free parameters. Taking advantage of the flexibility of the accelerator a delib- 
erate study of energy- and system-size dependence of the new phenomena is underway. The lively 
cooperation of theorists and experimentalists (strengthened by some new initiatives) will certainly 
lead to our much deeper understanding of the strongly coupled quark-gluon plasma. 
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